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Abstract 

During the last 35 years three giant flares were observed from so-called Soft 
Gamma Repeaters (SGR’s). They are assumed to be associated with star-quakes 
of pulsars accelerating electrons and, possibly, protons to high energy in the huge 
magnetic fields as inferred from the observations. Because of this and the observa¬ 
tion of non-thermal emission it has been speculated that they may be cosmic ray 
accelerators producing gamma-rays up to TeV energies. Neutrino telescopes, such 
as AMANDA and the ANTARES now under construction, could be used as TeV-7 
detectors for very short emissions by measuring underground muons produced in 7 
showers. We estimate signal and background rates for TeV photons from SGR giant 
flares in AMANDA, and we provide an estimate of the gamma shower events that 
Milagro could detect. 

Moreover, we consider that, if hadrons are accelerated in these sources, high 
energy neutrinos would be produced together with photons. These may be detected 
in neutrino telescopes using neutrino-induced cascades and upgoing muons. 

We argue that the Antarctic Muon and Neutrino Detector Array AMANDA may 
have observed the December 27, 2004 giant flare from the soft gamma-ray repeater 
SGR 1806-20 if the non-thermal component of the spectrum extends to TeV energies 
(at present the actual data is subject to blind analysis). Rates should be scaled by 
about two orders of magnitude in km^ detectors, such as IceGube, making SGR 
flares sources of primary interest. 


* Gorresponding author. 

Email addresses: halzen@pheno.physics.wisc.edu (Francis Halzen), 
hagar. landsmanSicecube. wise.edu (Hagar Landsman), 
tmontaruli@icecube.wisc.edu (Teresa Montaruli). 


Preprint submitted to Elsevier Science 


2 February 2008 





1 Introduction 


Soft gamma repeaters are X-ray pulsars that sporadically emit energies typ¬ 
ically of the order of ergs in X-rays and soft 7-rays, with the dis¬ 

tance D in units of 10 kpc, in short bursts lasting hundreds of milliseconds 
[1,2], Their steady emission is periodic with typical periods of 5-10 seconds 
and luminosities between 10^® -7 10^® ergs/s. These objects are thought to be 
magnetically-powered spinning neutron stars also known as magnetars. 

The first SGR was observed on Jan. 7, 1979 when a burst of soft gamma-rays, 
lasting about 0.25 s, was detected by the Venera space-craft from SGR 1806- 
20 [1]. To date five SGR’s have been identihed. The giant bursts are thought 
to result from the rearrangement of the huge magnetic held inferred to be of 
the order of R ~ 10^^ Gauss from their spin down rates [3]. Star-quakes are 
thought to fracture the rigid crust and cause these outbursts. The hares result 
from the formation and dissipation of strong localized currents resulting from 
magnetic held rearrangements associated with the quakes. 

SGR’s have been classihed as separate phenomena from gamma-ray bursts 
(GRB’s) and anomalous X-ray pulsars (AXP’s) [4] despite similarities. SGR 
bursts are less frequent than GRB’s and, unlike GRB’s, they do repeat and 
show clear periodicity. It is not excluded however that at least some of the 
short GRB’s are extra-galactic SGR’s and, vice-versa, that some of the intense 
hares of SGR’s are galactic GRB’s. SGR’s are very similar to the 8 observed 
AXP’s, the main diherence being the lack of bursts for the latter. Recently, 
though, SGR-like bursts were observed from at least one, and possibly two 
AXP’s [4]. Hence SGR’s and AXP’s might belong to the same class of neutron 
stars. 

Like GRB’s, SGR’s also produce non-thermal radiation during the short bursts. 
Because of this and the large luminosities produced, it has been speculated 
that these objects accelerate cosmic rays to high energy. The sudden changes 
of the large magnetic helds would accelerate protons or nuclei that produce 
neutral and charged pions in interactions with thermal radiation. These would 
be the parents of gamma-rays and neutrinos reaching high energies. 

We here propose that these high energy gamma-rays can be detected by un¬ 
derground nentrino telescopes and by experiments, snch as Milagro [5] and 
the Tibet array [6]. (The premier instruments for TeV gamma-ray detection, 
air Gherenkov telescopes, must point at an SGR in order to detect it and this 
is nnlikely for such a serendipitous event.) Undergronnd detectors have the 
capability to detect high energy mnons produced in air showers initiated by 
TeV-energy gamma-rays at the top of the atmosphere. Though gamma-ray 
showers are muon-poor, the large nnmber of secondary photons produced will 
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occasionally produce pions rather than electron-positron pairs. These decay 
into muons that can reach large underground detectors such as Baikal [7], 
AMANDA/IceCube [8,9] and future Mediterranean detectors [10,11,12,13]. 

TeV gamma-rays could be produced in electromagnetic process and also in 
hadronic interactions on the thermal radiation, and in this case they would 
be accompanied by a flux of neutrinos. A model for neutrino production from 
magnetars in their steady phase of periodic emission has been proposed [14]. 
Here we concentrate instead on the possibility that short bursts of neutrinos 
are produced during the rare violent giant flares. Because of the negligible 
background during the short time and in the direction of the SGR, the instru¬ 
ments have a much larger sensitivity for such events. 

This paper is organized as follows: in Sec. 2 we review the photon observa¬ 
tions of the giant flares. In Sec. 3 we argue that they should produce detectable 
events in earth-based experiments provided their spectra extend to TeV ener¬ 
gies not accessible to satellites. For illustration we calculate the rates of muons 
from gamma showers on Milagro surface [5] and the number of muons from 
gamma showers reaching the AMANDA detector [8]. In Sec. 4 we estimate 
the number of upward-going muons and showers produced by neutrinos in 
the same detector. Similar rates are expected in neutrino telescopes of similar 
depth such as ANTARES [10]. 


2 Giant bursts from Soft Gamma Repeaters: Observations 

In recent years SGR giant bursts have been observed from 3 objects: SGR 
0526-66 on May 5, 1979, SGR 1900+14 on August 27, 1998 and SGR 1806- 
20 on December 27, 2004 [1]. Their UT times and equatorial coordinates are 
given in Tab. 1. 


Source 

RA 

(J2000) 

Declination 

(J2000) 

Distance 

(kpc) 

UT 

SGR 0526-66 

05'*26™00.89^ 

-66°04'36.3" 

50 

May 5, 79 15:51 

SGR 1900+14 

19'*07™14.24^ 

+09° 19'20. l" 

14.5 

Aug. 27, 98 10:22:15.7 

SGR 1806-20 

18'*08™39.32* 

-20°24'39.5" 

15.1 

Dec. 27, 04 21:30:26.65 


Table 1 


Coordinates and UT times of the giant flares [1]. UT time for SGR 1900+14 is from 
[15] and for SGR 1806-20 from [16]. Positions (right ascension RA and declination) 
and distances are from [17] for SGR 1900+14 and from [18] and [19] for SGR 1806- 
20 , respectively. 

Giant flares are distinguished from common SGR bursts by their extreme 
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energies ergs/s) emitted during the first second followed by subse¬ 

quent emission lasting several minutes showing pulsations associated with the 
neutron star. While the total energy emitted is 6-7 orders of magnitude less 
intense than for GRB’s, their distance is typically a factor of 10^ smaller re¬ 
sulting in a flux increased by a factor 10^° relative to cosmological GRB’s. 
The opportunity for a detection is obvious. 

The bulk emission of SGR’s can be routinely explained by Optically Thin 
Thermal Bremsstrahlung (OTTB). We focus on the non-thermal component 
observed in giant flares with a hard power law (PL) spectrum, ^ oc E~^-^ 
[1,15]. This is harder than the Band spectrum observed for a typical GRB 
[20]. Because of the energetics and the non-thermal spectra giant flares should 
be prime suspects for accelerating cosmic rays and for producing associated 
fluxes of photons and neutrinos of TeV-energy, possibly higher. 

An interesting feature is that both SGR 1900-1-14 and SGR 1806-20 are as¬ 
sociated with clusters of massive stars [1,17,19], while the other giant flare 
emitter, SGR 0526-66, is associated with the supernova remnant N49 in the 
Large Magellanic Gloud. It has been suggested that association of massive 
stars, such as Gygnus OB2, can host highly magnetized pulsars that acceler¬ 
ate heavy nuclei up to high energies [21]. TeV y’s and z/’s can be produced 
by cosmic rays by interacting with the innermost parts of the winds of mas¬ 
sive O and B stars [22]. Moreover, relativistic nuclei photo-disintegrate in the 
intense radiation helds producing neutrinos from the decay of neutrons [23]. 
For a review of galactic neutron stars in high density environments: see [24]. 
In a recent paper [25] it has also been suggested that the December 27, 2004 
flare has similar features to a GRB and can be similarly explained as a rel- 
ativistically expanding hreball with baryon loading. Proton acceleration will 
inevitably lead to neutrino emissions. 

The December 27, 2004 flare exceed previous flares in intensity by over one 
order of magnitude. It was initially detected by INTEGRAL [16], then by 
Konus-Wind [26] and later by RHESSI [27]. While the germanium detectors 
operating from 3 keV to 15 MeV were saturated, RHESSI measured the flux 
using the particle detector with a time resolution of 0.125 s in two energy 
channels with thresholds of 65 keV and 650 keV. These data indicate signifi¬ 
cant emission above 650 keV for ~ 0.25 s during the giant flare. A conservative 
lower limit on the total fluence of the primary giant peak is 0.1 ergs/cm^, more 
than one order of magnitude larger than reported for SGR 1900-1-14 (5.5 • 10“^ 
ergs/cm^ in the first ~ 0.35 s [28]). Assuming a distance of ~ 15 kpc for 
SGR 1806-20 [19] and isotropic emission, a lower limit on the total hard X- 
ray/gamma-ray energy released is 8 • 10^® ergs during the giant flare. 

So far no spectral analysis has been reported for the 27 Dec., 2004 giant flare 
during which the AMANDA neutrino telescope was taking data. In Sec. 2.1 
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we study instead the features of the giant flare SGR 1900+14 on August 27, 
1998 during which the AMANDA detector had been partially completed and 
and was taking data with 10 strings; we refer to it as AMANDA B-10. 


2.1 The features of the SGR 1900+14 giant flare 


SGR 1900+14 was first observed in Jan. 1979 in the Aquila constellation. 
Thought initially to be located near SNR G42.8+0.6 at ~ 5 kpc, it is now 
believed to be associated with a massive star cluster at 14.5 kpc [17]. The 
position of the persistent emitter is given in Tab. 1 [17]. Its period of steady 
X-ray emission is 5.16 s with a luminosity of the order of 2 — 3.5 • 10^^ ergs/s. 
The quiescent energy spectrum can be htted by a Black Body (BB) function 
with a temperature of 0.43 keV and a Power Law (PL) shape with a negative 
spectral index in the range —1.0-j—2.5 [1]. Various flares with peak luminosi¬ 
ties between 10^® + 10^^ ergs/s and durations of the order of 0.1 s have been 
observed. SGR 1900+14 emitted a giant flare on Aug. 27, 1998 [15,28,29,30]. 
Two relatively high fluence events occurred on Aug. 29, 1998 [31] and on Apr. 
28, 2001 [32]. For reasons already explained in Sec. 2, we focus on the giant 
flare whose spectral properties have been well stndied, inclnding those of the 
hard non- thermal component [15,30]. 

Most information on the giant flare comes from Beppo-SAX [15,30], Konus- 
Wind [28] and Ulysses [29]. All snffer from satnration problems dnring the first 
second of the bnrst. Unfortnnately, there are no clear spectral measnrements 
covering the first ~ 0.35 s above 250 keV. Konus-Wind reported spectra below 
250 keV for the first second of the emission; see Fig. 1 [28]. Also shown is the 
connting rate showing the satnration of the detector. Konus-Wind spectra 
can be described by a OTTB spectrnm typical of electromagnetic processes 
because it was not sensitive to the high energy component. 

The energy and flnence measurements are also complicated by pile-up and 
dead time problems. Mazets et al. [28] hnd a lower limit to the flnence in the 
hrst ~ 0.35 s of 5.5 • 10“^ ergs/cm^ for > 15 keV. For a sonrce at 14.5 kpc 
and isotropic emission this corresponds to a lower limit on the total energy 
emitted by the source of 1.4 ■ 10^"^ ergs during the peak of the burst. A similar 
amount of energy was emitted in the long tail with an exponential decay of 
~ 90 s that followed the burst and lasted several hnndreds of seconds. 

Similar conclusions can be inferred from the Ulysses data [29] that were even 
more affected by satnration because of its sensitivity in the interval 25-150 keV. 

We next concentrate on the Beppo-SAX data with energy np to 700 keV. 
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Fig. 1. On the left: Spectra measured by Konus-Wind in 4 time intervals of the 
initial emission. On the right: Counting rate during the first second of the burst 
for Ery > 15 keV and parameter kT of the OTTB law that fits the 

energy spectra on the left as a function of time. Figures from [28]. 

2.2 Beppo-SAX observations of the Aug. 21, 1998 giant flare 


Beppo-SAX data on SGR 1900+14 were first published in [15] and updated 
in [30]. Data from the Beppo-SAX Gamma-Ray Burst Monitor (GRBM) were 
also affected by saturation and only a lower limit on the fluence is reported for 
the hrst 68 s of the flare. The experiment could not determine the event rise¬ 
time for timescales shorter than 1 s. Fig. 2 shows the counting rate as well as 
the energy spectra corresponding to 3 different time intervals in the 70-650 keV 
interval. Most importantly, the experiment provides information during the 
hrst 68 s where the data exhibit a very hard non-thermal component. The 
results were updated in [30] for a better understanding of the response of the 
GRBM. The energy range was extended to 40-700 keV. 

Possible hts to the data are presented in Fig. 3. The published best ht is 
a sum of 3 functions given in Tab. 2. We have instead chosen to adopt the 
previous ht in Feroci et ai, also shown in Fig. 3, because it gives a better 
description of the high energy component to which we are mostly interested 
in this work. The data are ht with OTTB + PL, while in the more recent 
paper a thermal bremsstrahlung spectrum (based on [33])+ BKNPL + PL 
(with positive spectral index, probably an artifact of the large errors in the 
interval htted) was adopted. The of fhe OTTB+PL is worse because the 
low energy part of the spectrum is poorly htted by the 2 functions compared 
to the BREM+BKNPL+PL ht. We also note that in the updated analysis the 
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systematic error was reduced from 10% to 2% in [15] in the low energy range 
and the minimum considered energy was decreased from 70 keV to 40 keV. 

The total fluence corresponding to the 3 intervals is: > 6.4 for interval A ^, 
2.1 for interval B and 0.49 for interval C, in units of 10“^ erg/cm^. Assuming 
an isotropic emission, the energy emitted by a source at 15.1 kpc is: > 1.7, 
0 .6, 0.1 in units of 10^^ erg. As already observed, the value in the hrst interval 
is smaller than the one reported by Konus-Wind for the first 0.35 s [28]. 

For further calculations we will adopt a differential energy spectrum for E > 
250 keV of ^ oc . We will study variations of the PL spectral index to 
allow for the uncertainties; see Fig. 4. Moreover, since typical Fermi accelera¬ 
tion processes are characterized by an E~‘^ spectrum, and given the snggestion 
that SGR 1806-20 may be a soft GRB [25], we also consider this assnmption 
for the spectral shape. Our different assnmptions are summarized in Tab. 3. 
Also given is the fraction of the fluence in the PL component relative to the 
total flnence measured by Beppo-SAX. From Fig. 4 it is clear that the different 
fnnctions adeqnately £t the data that were initially fitted as (16- (i7/keV)“^'^'^ 
keV“^ cm“^ s“^). For each spectral shape, the normalization was chosen so 
that the flnence in the interval 200-500 keV is the same as for the E~^'^'^ fit 
in [15]. 


1.1 • lO^E 40.3 + 

3.8 • 10-4^1-79 (E < 71keV) 

2.2 • 10®E-3-48 (E > 7lkeV) 


8.0 • IO-ioE+214 

BREMS 

BKNPL 


PL 

192.8E-%- 

■E/55 _|_ 

16.0E-i'^^ 

OTTB 


PL 


Table 2 

Functions used by Guidorzi et al. (BREMS-|-BKNPL-|-PL) [30,34] and by Feroci et 
al. (OTTB-I-PL) [15,34] to fit the spectrum for the first 68 s of the flare, respectively 
in the energy ranges 70-650 keV and 40-700 keV (see Fig. 3). Energy is in keV. 


3 TeV-Gamma Showers and Induced Underground Muons 


High energy gamma-rays initiate electromagnetic showers in the atmosphere 
via the dominant processes of bremsstrahlnng and pair prodnction. With the 
development of the shower a large population of photons emerges that have a 
probability to prodnce a pion rather than an electron pair. Snch pions decay 
into muons that, with sufficient energy, can penetrate to the depth of an 

^ Notice that this number is obtained as the integral of the functions in Tab. 2 
multiplied by the duration of the burst At = 68 s. 
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Fig. 2. On the left: Beppo-SAX counting rate during the SGR 1900+14 flare: 
the 3 intervals where the spectral analysis is done are 0-67 s (interval A), 68-195 s 
(interval B) and 196-323 s (interval C). On the right: The spectra for the 3 time 
intervals in the energy range of 40-700 keV. The best fit functions are shown too 
(figures from [15]). 


GF98 A GFSS A: 0rrB(kT=55 keV) -h PL(alDha=1.47) 




Fig. 3. On the left: Best fit to the spectrum measured by Beppo-SAX during 
the first 68 s of the flare (referred to as GF98 A in the plots) in the 40-700 keV 
range. On the right: fit using the Feroci et al. OTTB+PL function [15,34]. The 
corresponding fit functions are given in Tab. 2. Courtesy of C. Guidorzi [34]. 

underground detector. The fact that the process is rare is compensated by the 
large area of a detector such as AMANDA [35,36,37]. The showers can also be 
directly observed by ground-based wide held of view gamma detectors such as 
Milagro [35]. 


The high energy muons can be detected by deep under-ice/water neutrino tele- 
































Fig. 4. On the left: Upper curve consists of the OTTB with kT = 55 keV + the 4 
different PL components. Units are the same as in Beppo-SAK paper, keV“^ cm“^ 
[30]. This plot shows that the different PL components do not alter much the total 
spectrum at low energy. As a matter of fact all different PL’s added to the OTTB 
component overlap (see upper curves). On the right: Detail of the PL functions in 
the units used in this paper of keV“^ cm“^ s“^ as photon differential energy spectra 
for further estimates. 


Energy Spectrum 

Fraction of measured fluence in 68 s 

(keV~^ s“^ cm“^) 


0.2 X (E/keV)-^-'^^ 

12.9% 

16.0 X (E/keV)-^'^'^ 

12 .2% 

140.0 X (E/keV)-^-^^ 

14.8% 

329.1 X (E/A:eU)-2-00 

16.0% 


Table 3 


Differential photon energy spectra assumed in this paper. The fraction of the total 
fluence 6.4 • 10“'^ ergs measured by Beppo-SAX. between 40-700 keV that the PL 
components would account for is also reported. Numbers in bold are for the function 
used in [15]. 

scopes above the muon threshold determined by the detector depth, geometry 
and reconstruction capabilities, as observed in [35,36,37]. In order to derive 
the number of muons from a given photon spectrum we follow the analytical 
estimate in Ref. [38], that was checked against the Monte Carlo calculation 
in Ref [39]. We repeat here the main steps of that calculation. The muon 
spectrum from a differential photon spectrum ro(T^) is given by: 


dN^ 

dE 


ToiE) 


1 — 


A. 


EcosO \3{1 — r)/ - 


1 + In ■ 


A, 


( 1 ) 


where energies are in TeV, 6 is the zenith angle, r = rn^^/rn^ is the ratio of 
the squared muon and pion masses, e^r = 0.115 TeV is the energy above which 
pion interaction dominates over decay. The pion flux at the atmospheric depth 
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t due to the photon flux Fq can be factorized as n{E,t) = TQ{E)n 2 (t), with 
TTsW = r\j various terms are respectively: the 

z moment z^t^ = Jq that depends on the spectral index a 

of the photon flux ^ oc the effective pion interaction length in the 

atmosphere A^; the interaction length associated with tt photo-production 
\n- Following the discussion in Ref. [38] we assume: 


Att 1 Acr^/V" 

~ z -~-X < nx 

A^N ^ttA 


( 2 ) 


with z^T, = I and z^^t, = | (this is the z moment for pion regeneration in 
the atmosphere). These values are suitable for hard spectra with a around 1. 
We assume a constant pion photo-prodnction cross-section ~ 0.1 mb and 
an average mass nnmber A = 14.5 in the atmosphere. This is a reasonable 
assnmbtion in the photon energy range ~ 0.1 — 10^ TeV [40,39] ^ . Moreover, 
we assnme ~ 198 mb [38]. The logarithmic term in eq. 1 depends on tmax, 
the depth where the photon energy in the cascade has become too low to 


prodnce muons of energy E. It is assnmed tmax = A/j/n 
Xr is the radiation length in the atmosphere. E/ < Xy. 


c 


i> 


where 


> is the 7 energy 

reqnired to produce a muon of energy E and we assume < Xy^^ >~ 0.25. 

480 mb. 


Moreover, in eq. 1 we have 


— Aa 


= (1 


Z-n 


with aji 


Att Att '' ^ O'R 

The above muon flnx is valid for energies large enough that muons do not 
decay, e.g. in our case, in which muons must have enough energy to survive 
underground. 


The photon spectra motivated in Sec. 2.2 are expressed in the general form 
Tq{E) = ^ = ^ ■ 10“^^ with the nnits and parameters given in Tab. 3. 
We assume that the spectra we get from SGR 1900-1-14 data are similar for 
SGR 1806-20, a burst with similar features and duration. Because SGR 1806- 
20 is more energetic by over one order of magnitude, this represents another 
conservative estimate. Moreover, we recall that the measured values are lower 
limits on the fluence for both bursts. We assume that the photon spectrum 
extends to Ey^max and in Tab. 4 we estimate the number of events for two 
values, 200 and 500 TeV, and show the number of signal events as a function 
of Ey^max in Fig. 5. Photons of energies in excess of ~ 500 TeV shonld be 
attennated over galactic distances dne to interaction on backgronnd photons. 
The minimnm mnon energy E^^min that muons should have to cross the ice 
overburden depends on the inclination of the sonrce respect to the detector. 
For AMANDA, located at the Sonth Pole, a simple relation connects the zenith 
angle and the declination of the source, that is 6* = 5 -|- 90° = 70°. Being the 

^ The resulting cross section from the FLUKA transport and interaction code [41] 
is in agreemet with this value of the yV cross section, as well as data in Refs. [42] 
at GeV energies. 
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detector at a depth of 1500 — 2000 m, a muon induced by gammas from the 
source will have to transverse about 5.1 km of ice overburden before reaching 
the detector. We calculate the survival probability of muons with energy larger 
than 10 GeV at AMANDA depths using the MUM code for high energy muon 
energy losses [43]. In this calculation we account for the muon spectrum and 
the parent gamma spectrum in Tab. 3. Hence we derive the number of muons 
at AMANDA depth and for the given source zenith as: 

/ dN 

-j^Psurv{depth,Efj)dEfj, (3) 

F. 


assuming that > Egamma,max ~ 0.25E^^max, for depth equal 

to 5.1 km. 


Muons in AMANDA in 0.3 s 



Fig. 5. Number of muons induced by gamma showers in 0.3 s for an almost horizontal 
area of 30000 m^ for a detector like AMANDA-II at a depth between 1.5-2 km for 
the 4 different photon sectra indicated and as a function of the maximum photon 
energy. 


The main background source to the muons produced by gammas emitted in 
the flare are atmospheric muons. In the horizontal direction the atmospheric 
muon rate is reduced by roughly 2 orders of magnitude with respect to the 
vertical direction. Therefore the background is about 1 Hz. Furthermore, since 
the bulk of the energy is emitted in only 0.3 s from a point-like source, the 
background is limited to the time interval and the direction of the burst. 
Conservatively considering a search bin of 8°, well above the resolution of the 
detector and comparable to the size of the declination band [44] chosen by the 
experiment after a signal to noise optimization, the estimated background rate 
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from atmospheric muons is at the level of 7r(8*7r/180)^/(27r) xO.3x1 = 3-10 ^ 
events in 0.3 s. Hence it is negligible. 

Where the signal is concerned, we tabulate in Tab. 4 the expected number 
of muons from the burst in 0.3 s for E^^max = 200 TeV and 500 TeV. We as¬ 
sumed a near horizon effective area close to the geometrical area of AMANDA 
of about 30000 m^. The effective area includes cuts on the quality of the muon 
reconstruction, and was estimated for harder muon spectra than atmospheric 
muons for which an average value is 22000 m^. In Fig. 5 we show how the 
results vary for different E^^max values and for the 4 spectra we considered. 
Similar event rates might be observed by the future ANTARES detector, lo¬ 
cated at 2050 — 2400 m in the sea, that would be able to observe muons 
and showers induced by neutrinos from another flare from the source, SGR 
1900-1-14, for more than 40% of a day. It will also be able to observe neutrino 
induced upward going muons and showers from SGR 1806-20. 

As expected, the event rate is a strong fnnction of the spectral index and of 
the normalization. In the absence of a possible detection, it is straightforward 
for AMANDA to constrain the high energy flux associated with the possibility 
that SGR’s are high energy accelerators, possibly similar to GRB’s. 


F, 

Spectral 

Index 

g Events 

AMANDA 

Atmospheric 

//’s in AMANDA 

7 showers 

in Milagro 

6.47 • 10^3 

-0.73 

0.1 - 1 • 10 ® 

3• 10 -® 

1 • 10 ^^ 

9.42 10® 

-1.47 

0.8-6 

3 • 10"® 

8 • lO'^ 

3.13 • 10® 

-1.85 

0.1 - 1 • 10“2 

4• 10 “® 

79 

3.29 • 10® 

- 2.00 

0.1 - 1 • 10“3 

2 • 10 "^ 

5 


Table 4 


The normalization of the 7 flux and the spectral index for the giant flare. The 
normalization units are such that the differential flux ^ ■ [FjTeV)°^ • 

10“^^ is in cm“^ s“^ TeV“^. The values of the Beppo-SAX fit are in bold. The 
corresponding number of muons in a detector like AMANDA-II of 30000 m^ in 0.3 
s is given. Two numbers are given in the third column; the smaller is obtained 
assuming a maximum photon energy of 200 TeV and the larger assuming 500 TeV. 
Also the atmospheric muon background in a search solid angular bin of 8 ° and in 
0.3 s is indicated. The gamma shower events in Milagro in 0.3 s are also shown for 
a giant flare in its field of view (see the text). 


3.1 Estimate of photon induced showers in Milagro 


SGR 1806-20 was located above the Milagro horizon at the time of the Dec. 27 
giant flare at a zenith angle of about 68°. The source might have been visible 
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for Milagro though we understand data at such zenith angle are not routinely 
analyzed by the experiment whose range covers 0 — 60 degrees. Nonetheless, 
we consider the effective area presented in [45,37] as function of the 7 energy: 
A{Er^) = /TeVY in cm^, where /3 = 1.39 for E^ > 1 TeV and {3 = 2.35 

for E^ G [0.3,1] TeV, and we use it for estimating the number of shower events 
in Milagro from an SGR burst eventually happening in a favorable sky region. 
Using the various gamma spectra we estimate event rates in Tab. 4 using the 
following formula: 


N. 


showers 


T 


-E^7,maa: = 100TeV 



-^7',min=0.3TeV 


dN^ 

dE^ 


X A{E^)dE^. 


( 4 ) 


For some of the spectral indexes Milagro may have seen a large signal during 
the first instants of the flare. 


4 Neutrino initiated showers and npward-going muons from giant 
SGR flares 


The photon fluxes previously discussed are accompanied by neutrinos because 
charged pions are produced at the source along with the neutral pions that 
decay into photons. Assuming that photons and neutrinos have this common 
origin, the neutrino spectrum can be derived from the photon spectrum by 
energy conservation along with simple kinematics [36,47]. Photons may cas¬ 
cade in the source to emerge with a steeper spectrum and they may even be 
absorbed, therefore only a lower limit on the neutrino flux can be derived 
using this technique. 

It is likely that the target for neutrino production are the thermal photons 
surrounding the pulsar. We therefore assume p — 7 interactions with the ambi¬ 
ent radiation around the neutron star. There may also be a contribution from 
neutron decay from photo-disintegration of nuclei. The relation between the 
two spectra is: 

E-/,max 

rp I rp 

-^7,min ^u,rmn 


where K = 1 for pp interactions and 4 for p — 7 . In pp interactions charged and 
neutral pions are produced in equal numbers yielding 2 7 ’s carrying 1/2 of the 
neutral pion energy, that is 1/3 of the proton energy, and so E^ ~ Ep/Q. If 
muons also decay, there are 2 muon neutrinos per 7 and so Ejj ~ Ep/12. For 
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p — 7 interactions, neutrinos and photons are produced from the A resonance 
decay: p + 7 ^ A ^ and p + 7 ^ A —>• pvr^ ^ 

P 77 . Hence E^rnin is determined by the threshold of the process for protons 
[47]: 


Ep = T 


{2mp + 


2mp 


2 m„ 


r X 1.23GeV 


( 6 ) 


where we will consider a Lorenz factor L ~ 1. The minimum neutrino and pho¬ 
ton energy are about 5% and 10%, respectively, of the proton energy {E^^min = 
7 Xp — y'y 7 Ep^yyiiYi/^ uud E^prnin ^ —>7 ^ EppYiiyi j22j^ where 7 Xp — 0.2 is 
the average fraction of proton energy transferred to pions). The maximum pro¬ 
ton energy, and hence the maximum gamma and neutrino energies, depend on 
the kinematics in the beam that results from the burst. We will assume values 
a bit larger than what is calculated in Ref. [14]: here Ep^^nax ~ 100 — 300 TeV 
[14], where protons are accelerated in the potential drop of the magnetar with 
surface magnetic held of ~ 10 ^® Gauss. As a matter of fact, we are considering 
giant hares and not the steady emission. Hence, we will assume here a maxi¬ 
mum photon energy of 200 TeV and a maximum neutrino energy of 100 TeV. 
Also we assume that the photon spectra have the same spectral index, a very 
conservative assumption [37]. Using the 7 spectra in Tab. 3, we obtain the 
neutrino huxes in Tab. 5. 


Neutrino Spectrum 

(cm“^ GeV“^) 

Cascade events (s 

l^e Vp Vr 

Upward-going /x’s 

(s-i) 

5.90(E/GeV)-°'^3 

1 • 10^ 

8 • 10® 

1 • 10^ 

2 • 105 

8.74 • 10-3(7;/GeU)-i-^^ 

0.4 

0.2 

0.3 

0.8 

3.09 • 10-^(A/GeU)-i-85 

9 • 10-5 

3 • 10-5 

1 

0 

r —1 

0 

5 • 10-4 

8.23 • 

3 • 10-5 

1•10-5 

2 • 10-5 

3 • 10-5 


Table 5 

Neutrino spectra obtained using the fluxes in Tab. 3 and eq. 5 and shower and 
upgoing-muon event rates in a detector like AMANDA-II per second for a source 


at (5 = ±20°. 

We use these huxes to derive event rates for two diherent neutrino signatures: 
neutrino induced showers (cascade events), that have a dvr coverage and can 
be observed even if the source is in the upper hemisphere of the detector and 
neutrino induced upward-going muons that can be observed only if the source 
is in the lower hemisphere; this is the case for SGR 1806-20 in the held of 
view of ANTARES and for SGR 1900±14 observable by AMANDA-H. We 
estimate event rates for AMANDA neutrino ehective areas, though similar 
numbers apply for the present ANTARES detector design. Neutrino ehective 
areas as a function of the energy for the declination of a source can be directly 
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convoluted with the neutrino spectra to provide event rates. To estimate the 
number of cascade events we have considered the AMANDA effective areas 
for neutrinos of all flavors presented in [48]. Because of oscillations, a cosmic 
neutrino beam reaches the Earth with equal fluxes for all flavors. For the 
near horizontal direction of the sources under consideration the absorption of 
neutrinos is almost negligible. To estimate upward-going neutrino event rates 
we used the effective area presented in [44]. Results are given in Tab. 5. 

The background for neutrinos from the flare is due to atmospheric neutrinos 
and muons. As already explained in Sec. 3, it can be rejected using time and 
directional constraints. We neglect the contribution to the neutrino signal from 
the flare due to neutrinos that accompany muons from gammas considered 
in Sec. 3 because they would give a much smaller contribution to neutrinos 
directly produced at the source. 


5 Conclusions 


Starting from measured spectra by Beppo-SAX for the SGR1900-I-14 giant 
flare, we have estimated the TeV gamma-ray and neutrino event rates ex¬ 
pected (and possibly observed) for the AMANDA detector from SGR 1806-20 
on Dec. 27, 2004. We showed that tight time and angular constraints can re¬ 
duce the atmospheric muon background to negligible levels. The detectors are 
essentially free of background because of the time constraints. 

A negative result will result in an interesting constraint on the possibility that 
SGR most intense flares are similar to small GRB’s in our Galaxy. 
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